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I.  INTRODUCTION 


.V 


It  is  well  known  that  when  an  electromagnetic  field,  coupling  two  levels 
of  an  atomic  or  molecular  system,  has  its  frequency  swept  across  resonance, 
population  reversal  can  occur  as  a  result  of  the  coherent  transient  phenomenon 
known  as  adiabatic  rapid  passage  (ARP)*^^“  For  this  phenomenon  to  occur,  the 

rj. 

frequency  sweep  rate  is  required  to  be  rapid  with  respect  to  the  relaxation 
processes  occuring  in  the  system,  and  at  the  same  time  slow  enough  so  that  the 
system  can  follow  the  frequency  sweep^idiabatically.  Though  ARP  was  first 
observed  in  nuclear  magnetic  resonance  (NMR),  the  phenomenon  is  quite  general 
and  has  been  demonstrated  in  electron  paramagnetic  resonance  (EPR)  experiments 
and  in  both  one-^*^  and  two-^  photon  optical  studies.  i>ihe  present  investiga¬ 
tions  treat  ARP  in  the  ]^Rb  hyperfine  transition,  where  optical  pumping  pro¬ 
duces  a  population  imbalance  in  the  ground-state  hyperfine  levels,  5  and 
where  the  Influence  of  a  microwave  frequency  swept  across  the  hyperfine  reso¬ 
nance  is  monitored  by  the  transmitted  intensity  of  the  optical  pumping  laser 
Though  the  hyperfine  transition  of  ®^Rb  is  interesting  in  its  own  right,  our  y 
prime  concern  was  to  obtain  a  better  understanding  of  ARP  processes  in 
general. 


The  conditions  on  the  frequency  sweep  rate  for  ARP  are  typically  written 


as 


«  ||_  u  -  «0)|  «  Ul2  (l) 

where  u>j  is  the  Rabl  frequency,  T ^  the  dephaslng  time,  and  (o>  -  u>q)  the  detun¬ 
ing  between  the  resonant  frequency  and  the  frequency  of  the  perturbing  field. 
Defining  B  as  the  rate  of  change  of  the  resonance  condition,  Eq.  (1)  can  be 
written  in  the  dimensionless  form 


where  the  terms  on  the  left-  and  right-hand  sides  of  Eq.  (2)  specify  the  sweep 
requirements  of  rapidity  and  adiabaticity,  respectively.  Whereas  population 
reversal  is  typically  stated  as  the  result  of  ARP  when  these  two  conditions 
are  satisfied,  in  all  real  systems,  with  finite  values  of  T£,  total  population 
reversal  can  only  occur  when  the  conditions  are  maximally  satisfied,  l.e., 
when  the  Rabi  frequency  is  infinite.  In  the  laboratory  this  can  never  be 

O 

accomplished;  instead,  "we  take  as  nature  gives  it  and  work  with  some 
finite  value  of  a>j.  This,  however,  raises  a  question:  How  should  the  sweep 
rate  0  be  chosen  in  order  to  maximize  the  population  reversal?  Analysis  of 
this  question  in  the  Rb  system  and  application  of  the  results  to  the  generic 
ARP  process  is  the  thrust  of  this  paper. 

To  aid  in  the  analysis  of  ARP  two  parameters  are  defined.  the 

rapidity  parameter,  and  n^,  the  adiabaticity  parameter,  are  equal  to  the  terms 
on  the  left-  and  right-hand  sides  of  Eq.  (2),  respectively,  which  maximize 
population  reversal  during  ARP.  To  facilitate  general  applicability  of  the 
current  results,  these  two  parameters  are  analyzed  in  terms  of  a  dimensionless 
quantity,  the  normalized  Rabi  frequency  u^T^,  and  from  Eq.  (2)  we  see  that 

n0^0  “  “lT2  ^ 

Essentially,  we  want  to  determine  the  behavior  of  £q  and  as  functions  of 
the  normalized  Rabi  frequency,  and  to  understand  on  a  conceptual  level  the 
physics  behind  this  behavior.  Though  Eq.  (2)  states  the  general  conditions 
for  ARP,  it  gives  no  clue  as  to  the  behavior  we  seek;  furthermore,  we  have 
found  no  detailed  discussion  of  this  behavior  in  the  literature. 

It  is  particularly  interesting  to  study  ARP  in  the  ground-state  hyperfine 
87 

levels  of  Rb,  as  this  is  not  a  two-level  system.  The  two  hyperfine  levels 
87 

of  Rb,  F  ■  1  and  F  ■  2,  are  composed  of  Zeeman  sublevels,  of  which  only  the 
nip  ■  0  levels  are  coupled  by  the  sweeping  radiation  field.  In  addition  to 
experimental  study,  theoretical  analyses  are  conducted  at  several  levels. 
First,  a  detailed  density-matrix  analysis  of  the  entire  8^Rb  hyperfine  system 
Is  performed.  The  system  is  then  reanalyzed  in  terms  of  a  strict  two-level 


model.  This  analysis  provides  insight  into  how  applicable  the  current  experi 
mental  results  are  to  other  two-level  ARP  processes.  Finally,  a  simplified 
analysis  of  the  two-level  system  yields  an  approximate  closed-form  result  for 
the  degree  of  population  reversal  during  ARP.  From  this  result  we  feel  phys¬ 
ical  Insight  into  the  ARP  process  is  gained. 


II.  EXPERIMENTAL  PROCEDURE 


A.  GENERAL  DESCRIPTION 

In  the  present  experiment  we  Investigated  the  conditions  for  ARP  in  the 
ground  state  0-0  hyperfine  transition  of  ®^Rb;  Fig.  1  shows  the  radiative 
interactions  of  interest.  A  single-mode  diode  laser  (Mitsubishi  ML-4101), 
tuned  to  one  of  the  Dj  hyperfine  resonance  lines  at  794.7  nm,  was  used  to 
optically  pump  atoms  from  the  5^S  y^  (F  ■  2)  Zeeman  manifold  into  the  5^S  y 
(F  -  1)  Zeeman  manifold.  The  Doppler-broadened  absorption  linewidth  (-500 
MHz)  was  greater  than  both  the  laser  linewidth  (~100  MHz)  and  the  Zeeman 
splitting  (<700  kHz).  Thus  optical  pumping  occurred  from  only  the  F  *  2 

Q 

hyperfine  state,  but  from  all  Zeeman  sublevels  of  that  state. 

The  transmitted  intensity  of  the  laser  through  the  Rb  vapor  was  used  to 
determine  the  density  of  atoms  in  the  F  •  2  level,  and  in  steady-state,  in  the 
absence  of  the  microwave  field,  this  intensity  is  a  maximum.  By  sweeping  a 
microwave  field  across  the  0-0  hyperfine  transition,  atoms  in  the  (F  -  1, 
mF  ■  0)  state  were  transferred  into  the  (2,0)  state,  and  the  resulting  change 
in  the  transmitted  laser  intensity  was  a  measure  of  the  number  of  atoms  trans¬ 
ferred.  For  a  fixed  microwave  field  strength,  the  sweep  rate  of  the  microwave 
frequency  was  adjusted  so  we  could  observe  the  largest  change  in  the  trans¬ 
mitted  light  intensity.  Thus,  we  were  able  to  measure  the  sweep  rate  that 
maximized  the  population  reversal  as  a  function  of  microwave  field'  strength 
or,  equivalently,  the  microwave  Rabi  frequency. 

B.  DETAILS  OF  THE  EXPERIMENTAL  ARRANGEMENT 

The  experimental  apparatus  is  shown  schematically  in  Fig.  2.  A  Corning 
7070  glass  absorption  cell  that  contained  an  excess  of  *^Rb  metal  and  N2  at  10 
Torr  was  situated  in  a  TEqH  cylindrical  microwave  cavity  (R  -  2.8  cm)  tuned 
to  the  ®^Rb  ground-state  hyperfine  transition  frequency,  6834  MHz.  The  N2  was 
present  in  order  to  quench  the  Rb  fluorescence,  and  to  act  as  a  buffer  to 
reduce  the  effect  of  collisions  with  the  cell  walls.  A  static  magnetic  field 
of  a  few  hundred  milligauss  was  applied  parallel  to  the  cavity  axis  in  order 
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Fig.  2.  Experimental  Arrangement  for  the  ARP  Experiments.  In  order 
to  measure  the  slow-passage  linewidth  in  the  limit  of  zero 
microwave  power,  the  preamp  was  replaced  by  a  lock-in 
amplifier,  and  the  microwaves  were  chopped  at  a  low 
frequency. 
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to  define  the  quantization  axis,  and  to  split  the  Zeeman  levels  so  that  only 
the  0-0  transition  was  Induced  by  the  microwaves.  The  cavity  and  cell  were 
thermostatically  controlled  to  ±0.18C  at  about  37°C,  and  were  surrounded  by 
three  layers  of  magnetic  shielding. 

The  diode  laser  emission  was  collimated  by  a  short-focal-length  lens  to  a 
diameter  of  -0.8  cm,  so  that  an  aperture  (~0.2  cm  diameter)  allowed  only  the 
central  portion  of  the  laser  beam  to  enter  the  absorption  cell.  Since  the 
radial  profile  of  the  laser  emission  is  well  approximated  by  a  Gaussian,  and 
since  the  vapor  was  optically  thin  ([Rb]ol  ■  0.56),*®  the  intensity  distribu¬ 
tion  in  the  cell  volume  was  expected  to  be  fairly  uniform.  These  precautions 
were  necessary  in  crder  to  make  the  effects  of  light-induced  inhomogeneous 
broadening  negligible. *  * 

The  microwave  frequency  sweep  was  generated  by  applying  a  voltage  ramp  to 
a  calibrated  voltage-controlled  crystal  oscillator  (VCXO),  whose  output  at 
-106  MHz  was  multiplied  to  the  87Rb  hyperfine  frequency  region  and  then  atten¬ 
uated  before  being  allowed  to  power  the  microwave  cavity.  The  voltage  ramp 
amplitude  and  duration  determined  the  frequency  sweep  rate  of  the  6834-MHz 
microwaves,  and  these  were  adjusted  to  produce  a  maximum  signal  amplitude  on 
the  oscilloscope;  a  typical  ARP  signal  is  shown  in  Fig.  3.  The  sweep  rate 
that  produces  the  maximum  ARP  signal  is  designated  8q.  In  order  to  ensure 
that  our  measured  corresponded  to  a  single  pass  through  resonance,  starting 
with  the  system  in  a  steady-state  condition,  the  ramp  voltage  to  the  VCXO  was 
supplied  by  a  "slave"  function  generator  whose  output  was  triggered  every 
2  sec  by  a  "master"  function  generator.  Since  our  experimentally  determined 
pumping  and  relaxation  times  were  on  the  order  of  milliseconds,  steady-state 
conditions  were  attained  in  the  time  between  each  passage. 

As  a  result  of  the  presence  of  the  buffer  gas,  the  atoms  are  considered 

as  essentially  "frozen"  in  place;**  therefore,  they  experience  a  microwave 
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Rabi  frequency  determined  by  their  position  in  the  TEqjj  mode:4 


MICROWAVE  FREQUENCY 


The  Absorption  Signal  versus  Microwave  Frequency,  as 
Measured  by  the  Transmitted  Light  Intensity.  The  symmetric 
line  shape  corresponds  to  a  alow  passage  through 
resonance.  The  asymmetric  line  shape,  five  times  larger 
than  actually  shown,  corresponds  to  a  rapid  passage  that 
maximised  the  population  reversal  for  an  average  Rabl 
frequency  in  the  cavity  of  ~2.3  kHz.  Note  that  the  maximum 
of  the  fast-passage  signal  occurs  when  the  microwaves  have 
swept  approximately  one  half-width  beyond  resonance. 


where  R  and  L  are,  respectively,  the  cavity  radius  and  length.  The  experiment 
is  arranged  so  that  the  laser  beam  travels  down  the  center  of  the  cavity,  and 
because  the  laser  beam  diameter  is  so  much  smaller  than  the  cavity  diameter, 
the  radial  variation  of  the  Rabi  frequency  for  those  atoms  probed  by  the  laser 
is  only  ~0.5%.  Therefore,  the  Rabi  frequency  may  be  accurately  approximated 
by 


ui^r.z)  *  o>j(0,z)  -  o>lp  sin(|^-)  (5) 

Since  all  the  atoms  in  the  laser  beam  contribute  to  the  signal  amplitude,  we 
somewhat  arbitrarily  associate  our  measured  optimum  sweep  rate  with  an  average 
of  the  Rabi  frequency  over  the  length  of  the  cavity: 

<“1>  '  7  “ip  <6> 

<u)j>  is  then  the  Rabi  frequency  used  in  determining  the  normalized  Rabi  fre¬ 
quency  <*>jT2  associated  with  the  values  of  t^q  and 

The  absolute  magnitude  of  the  Rabi  frequency  within  the  cavity  is  not 
obtained  from  the  rapid-passage  experiment.  Rather,  the  quantity  over  which 
experimental  control  is  maintained  is  the  microwave  power  entering  the  cavity. 
Through  the  use  of  variable  microwave  attenuators,  precise  microwave  power 
levels  are  allowed  into  the  cavity.  To  relate  these  power  levels  to  absolute 
Rabi  frequencies  requires  additional  measurements. 

C.  MEASUREMENT  OF  THE  RABI  FREQUENCY.  THE  DEPHASING  RATE. 

AND  THE  RATIO 

It  is  well  known  that  the  half  width  at  half  maximum  (HWHM),  of 

the  transition  line  shape  in  a  two-level  system  is  given  by 

Am1/2-  (<1/T2)2  +  (Tj/T^  w2]1/2  (7) 

where  Tj  and  T2  are,  respectively,  the  longitudinal  and  transverse  relaxation 
times.  For  the  case  of  the  0-0  hyperfine  transition  in  an  alkali,  the  expres¬ 
sion  for  can  he  Put  *n  the  same  form  as  Eq.  (7).  However,  because  of 
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the  presence  of  the  other  Zeeman  sublevels,  Tj  is  no  longer  the  standard 

longitudinal  relaxation  time,  but  becomes  a  complicated  function  of  the  photon 

13 

absorption  rate  B  and  a  phenomenological  relaxation  rate  The  dephasing 

rate  I/T2  is,  however,  related  to  B  and  y  in  a  fairly  simple  way: 


(8) 


For  the  transmission  line  shape  in  our  experiment,  the  Rabi  frequency 

appearing  in  Eq.  (7)  is  well  approximated  by  u,  ;  this  is  discussed  in  the 

IP 

Appendix.  Line  shapes  were  recorded  on  a  strip  chart  for  each  attenuation 
level  of  the  microwave  power.  Care  was  taken  to  ensure  that  these  line  shapes 
corresponded  to  a  slow  passage  through  resonance.  Thus,  once  the  dephasing 
rate  and  the  ratio  (T1/T2)  were  determined,  we  could  extract  from  our 
linewidth  measurements. 

As  implied  by  Eq.  (7),  the  dephasing  rate  is  given  by  the  transmission 
linewidth  in  the  limit  of  zero  Rabi  frequency.  To  perform  the  measurement  we 
replaced  the  preamp  with  a  lock-in  amplifier  and  chopped  the  microwaves  with  a 
diode  switch  placed  after  the  microwave  attenuator.  In  this  way  we  could 
measure  linewidths  at  very  low  microwave  power  levels  and  extrapolate  these 
measurements  to  zero  microwave  power.  The  results  of  these  measurements  are 
shown  in  Fig.  A,  and  the  solid  line  is  a  least-squares  fit  of  the  data  to  the 
form  of  Eq.  (7).  This  procedure  yielded  an  extrapolated  dephaslng  rate  of 
52  Hz. 

In  order  to  determine  Tj/T2  it  was  necessary  to  know  both  B  and  y  indi¬ 
vidually.  This  was  accomplished  by  using  our  known  value  of  the  dephasing 
rate,  and  a  measurement  of  the  fractional  population  (  in  the  F  •  2  ground- 
state  hyperfine  level  as  a  result  of  optical  pumping.  The  transmitted 
intensity  of  the  laser  beam  was  measured  for  the  laser  tuned  off  and  on 
resonance : 


£n[I(off )/I(on)l  -  c  [Rb]  oi 


(9) 
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•  Plot  of  the  Measured  Half-Width  Half  —Maximum  (HWHM)  versus 
the  Square  Root  of  the  Noraalized  Microwave  Power  Fed  into 
the  Cavity.  The  solid  line  is  a  least-squares  fit  of  the 
data,  yielding  an  Intercept  of  52  Hz. 


L- 

>  m 


and  yielded  the  result  x,  *  0.31.^  The  homogeneous  density-matrix  rate  equa¬ 
tions  for  were  then  solved  in  steady  state, ^  and  the  pumping  rate  and 

relaxation  rate  were  varied  so  as  to  give  calculated  results  consistent  with 
our  experimental  determination  of  I/T2  and  5.  In  this  way  we  obtained  values 
of  60  Hz  and  22  Hz  for  B  and  y,  respectively,  which  resulted  in  ^1/^2  ■  1.34. 
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III.  RESULTS 


Using  the  various  values  of  and  It  Is  a  straightforward  procedure 


to  evaluate  the  corresponding  rapidity  and  adiabatldty  parameters.  The 


relevant  experimental  data,  and  the  resulting  determination  of  £q  and  t}q,  are 


collected  In  Table  1.  Figure  5  shows  these  values  plotted  as  functions  of 
normalized  Rabi  frequency,  Essentially,  this  figure  shows,  as  a  func¬ 

tion  of  Rabl  frequency,  how  rapid  and  how  adiabatic  the  passage  through 
resonance  must  be  In  order  to  maximize  the  degree  of  population  reversal. 


The  most  striking  feature  of  the  figure  is  the  disparity  in  the  rate  of 
change  of  the  parameters  as  a  function  of  Rabi  frequency.  The  adlabatlcity 
parameter,  while  greater  than  unity  (implying  satisfaction  of  the  adiabatic 
condition),  increases  very  slowly.  However,  the  rapidity  parameter  decreases 
at  a  very  much  faster  rate.  This  is  not  at  all  obvious  from  the  conditions 
for  ARP,  and  in  fact  we  had  naively  expected  the  rapidity  and  adlabatlcity 
parameters  to  change  at  approximately  the  same  rate.  The  implication  of  these 
results  is  that  while  it  is  important  to  satisfy  the  adiabatic  condition  in 
ARP,  once  it  is  satisfied  it  is  more  important  to  pass  through  the  resonance 
as  rapidly  as  possible. 


Table  1.  Relevant  Experimental  Data 


Microwave 
Attenuation 
Level ,  db 


rad/s2 


1/2  * 
rad/s 


1.1  x  10° 

6.7  x  107 

•>  5  w  in? 


2.5  x  10 ' 
1.3  x  107 

6.7  x  106 
3.0  x  106 

1.0  x  106 

5.7  x  105 


2.7  x  10“ 
1.9  x  104 
1.2  x  104 
8.5  x  103 
6.0  x  103 

3.8  x  103 
2.7  x  103 

1.9  x  103 


2.3  x  10“ 

i  i  w  mA 


1.7  x  10“ 
1.0  x  104 
7.4  x  103 

5.2  x  103 

3.3  x  103 

2.3  x  103 

1.6  x  103 


.043 

1 

.053 

1 

.083 

1 

.12 

1 

.16 

1 

.23 

1 

.48 

2 

.58 
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Plots  of  the  Adlabaticlty  Parameter  ru  and  the  Rapidity 
Parameter  as  Functions  of  the  Normalised  Rabi  Fre-  * 
quency  w^T..  The  boxes  and  triangles  correspond  to  our 
■easuremenfs  of  the  adiabaticity  and  rapidity  parameters 
respectively.  The  solid  line  represents  the  theoretical’ 
Zlll*  ^t8lned  *lth  the  ful1  ^  <*«nsity  matrix,  the 

ttl  tL ^^rfrfr8POnd<l  t0  th!  apin‘I/2  numerical  solutions, 
and  the  dotted  line  corresponds  to  our  approximate  formula. 

ilTthe  tlxtfe8ent  ***  daU  °f  Nakaniahl  et  al*»  •«  explained 
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IV .  THEORY 


We  have  considered  the  ARP  problem  within  three  different  theoretical 
frameworks:  (1)  a  detailed  density-matrix  approach,  where  we  include  relax¬ 
ation  among  the  various  sublevels  and  numerically  solve  the  resulting  coupled 
differential  equations;  (2)  a  two-level  approximation  of  our  system,  where 
the  coupled  rate  equations  are  solved  numerically;  and  (3)  an  approximate 
analytical  solution  to  the  two-level  problem.  The  first  set  of  calculations 
allows  us  to  make  a  comparison  with  experiment  in  order  to  determine  the 
validity  of  both  our  experimental  results  and  the  theoretical  model.  By 
comparing  the  second  set  of  calculations  with  the  first,  we  can  determine  the 
degree  to  which  our  observations  are  generally  valid.  The  last  calculation, 
though  approximate,  is  meant  to  provide  physical  Insight  into  the  theory,  an 
Insight  not  readily  available  from  the  numerical  solutions. 

A.  THE  87Rb  DENSITY-MATRIX  SOLUTION 

The  density-matrix  equations  used  in  the  present  calculation  are  essen¬ 
tially  the  same  as  those  used  previously  by  Camparo  et  al.,11  except  that  we 

do  not  need  to  consider  the  effect  of  light  shifts.  The  model  assumes  that 

2 

all  Zeeman  sublevels  of  the  5  (F  ■  2)  state  have  equal  optical  absorption 

rates,  and  that  fluorescence  repopulates  the  two  ground-state  hyperfine  levels 
with  equal  probability  at  a  rate  much  faster  than  does  optical  excitation. 
Furthermore,  it  is  assumed  that  only  two  phenomenological  relaxation  rates  are 
necessary  to  describe  the  ground  state’s  return  to  its  unperturbed  state:  a 
longitudinal  relaxation  rate  y ^  for  the  diagonal  density-matrix  elements,  and 
a  transverse  relaxation  rate  for  the  off-diagonal  density-matrix  elements. 
Lastly,  since  the  microwaves  are  only  in  resonance  with  the  0-0  hyperfine 
transition,  there  are  only  two  off-diagonal  elements  that  need  to  be 
considered. 

As  discussed  in  Ref.  11,  the  symmetry  of  the  model  and  the  normalization 
of  the  density  matrix  reduces  the  complexity  of  the  problem,  so  that  we  need 
only  find  solutions  to  five  coupled  differential  equations:^ 


with 


p(2,2) 

P(2,0)  y 

p  =  p(  1 ,0)  R  ”  ■g^- 

In  [p(2,0; 1 ,0)] 
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The  Instantaneous  detuning  between  the  hyperflne  resonance  and  the  microwave 
field  is  A(t),  where 

A(t)  -  AQ  -  Bt 

and  Aq  Is  the  detuning  at  time  t  ■  0;  we  assume  a  linear  change  in  the 
microwave  frequency  with  time :  du/dt  ■  B  • 

These  differential  equations  are  solved  numerically  using  Hamming's  modi¬ 
fied  predictor-corrector  method  with  a  fourth-order  Runge-Kutta  starter.^  He 
assume  y^  -  ^ »  ancl  employ  the  values  of  B  and  y  obtained  from  the  experiment. 
Aq  is  chosen  large  enough  so  that  the  equation  solutions  at  t  ■  0  can  be  well 
approximated  by  their  off-resonance  steady-state  values.  To  determine  the 
value  of  the  sweep  rate  that  maximizes  the  population  reversal  for  a  fixed 
Rabl  frequency,  the  differential  equations  are  solved,  yielding  the  quantity 

y(tft)  -  4  p(2,2)  +  p(2,0).  (14) 

U  C0  C0 


1 

1 

1 

0 
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for  various  values  of  £.  y(t)  is  the  total  density  of  the  F  =»  2  hyperfine 
level,  and  tQ  is  the  time  when  this  reaches  a  maximum.  An  example  of  the  set 
of  discrete  data  points  generated  in  this  way  is  shown  in  Fig.  6.  A  spline 
fit  is  then  used  to  Interpolate  between  these  points  in  order  to  determine  the 
theoretical  value  of  8q  for  the  particular  Rabi  frequency.^  The  theoretical 
values  for  the  rapidity  and  adiabaticity  parameters  obtained  in  this  way  are 
plotted  as  a  solid  curve  in  Fig.  5. 

B.  THE  TWO-LEVEL  DENSITY  MATRIX 

The  physics  associated  with  this  model  is  essentially  the  same  as  that 
described  in  the  previous  section,  except  that  we  now  consider  two  states 
rather  than  eight.  Letting  the  subscripts  1  and  2  denote  the  lower  and  upper 
hyperfine  levels,  respectively,  we  consider  solutions  to  the  three  coupled 
rate  equations: 

Y 

a22  "  T"  "  (f  +  Yl)°22  +  “l1®  (o21)  (I5a) 

Re(®21  >  “  “  (f  +  y2^  **^21^  +  A(c>  1111  ^°21*  (15b) 

and 

u . 

Im  (o21)  -  -  (y  +  T2)  “  A(t)  Re(°2p  +  d  “  ^22^  (15c) 

where  we  have  made  use  of  the  normalization  condition:  o 22  +  o  ^  *  1.  The 
procedure  for  the  numerical  solution  for  (5q  is  the  same  as  before,  except  that 
we  use  a  Runge-Kutta  technique  to  solve  the  equations.^  The  calculated 
values  of  the  rapidity  and  adiabaticity  parameter  for  this  model,  which  is 
essentially  a  Bloch  model,  are  displayed  in  Fig.  5  as  a  dashed  curve. 

C.  ANALYTIC  APPROXIMATION 

He  now  consider  a  spln-1/2  atom  whose  wavefunction  is  a  linear  combina¬ 
tion  of  the  spin-up  and  spin-down  basis  wavef unctions,  x+  and  x_: 

ip(t)  ■  Cj(t)X+  +  C2(t)X_  (16) 
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LOG  [Sweep  Rate  0  (rad/sec^)] 


I 


An  Example  of  the  Discrete  Set  of  Data  Points  Generated  by 
our  Model.  The  open  circles  correspond  to  the  signal 
amplitude,  In  arbitrary  units,  calculated  for  a  particular 
choice  of  the  sweep  rate  6.  The  solid  line  is  a  spline  fit 
of  the  numerical  data  points,  allowing  us  to  determine 


We  take  the  zero  of  energy  for  the  atom  as  midway  between  the  two  states  (E+  - 
E_  -  hu>g)  and  subject  the  atom  to  an  oscillating  electromagnetic  field  whose 
Instantaneous  frequency  is  changing  at  a  rate  6.  If  we  ignore  the  counter¬ 
rotating  component  of  the  electromagnetic  field,  the  wavefunction  satisfies 
the  Schrodinger  equation 


i  h  $>(t)  -  H( t)ifi( t) 


(17) 


where  the  Hamiltonian  is 

g  £ 

(coq  -  iy)  i^expl-i^-  |-)t] 

o)1exp[i(co0-  |^-)t]  -u)0 

and  where  we  have  included  relaxation  of  the  spin-up  state  in  the  form  of  an 
Imaginary  energy  term,  -  .  We  will  assume  that  for  times  very  much  before 

the  passage  through  resonance  the  atom  is  in  the  spin-down  state  (l.e., 
tHm-«  | C2 ( t )  |  ^  -  1  ) .  After  the  passage  through  resonance  there  will  be  some 
nonzero  probability  of  finding  the  atom  in  the  spin-up  state,  but  as 
t  -►  +»  this  will  decay  to  zero  as  a  result  of  the  damping  term  in  the 
Hamiltonian. 
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Horwitz  has  investigated  a  particular  form  of  this  problem,  that  of  y 
equal  to  zero.  He  found  that  in  this  case  |Cj(t)|^,  a  measure  of  the  popula¬ 
tion  reversal,  may  be  expressed  as  a  linear  combination  of  two  confluent 
hypergeometric  functions.  As  t  becomes  large  the  reversal  is  complete  and  the 
asymptotic  form  of  the  solution  may  be  used,  yielding 


t1*®*  |Cj(t)|2  -  [1  -  exp(-  ita)j/2g)]  (19) 

If  y  is  nonzero  the  problem  is  more  complicated.  However,  when  the  sweep  rate 
is  more  rapid  than  the  relaxation  rate,  as  it  is  for  the  cases  we  wish  to 
consider,  Horowitz  states  that  Eq.  (19)  is  only  modified  by  the  Inclusion  of 
an  exponential  decay  term.  Thus,  in  the  presence  of  relaxation  we  expect 

lci( tm) 1 2  •  U  “  exp(-  nwj/2g)]  exp(-  yt^)  (20) 


with  t0  the  time  corresponding  to  maximum  population  reversal 


To  specify  tm,  use  is  made  of  our  experimental  observations.  The  maximum 
of  the  ARP  signal  occurs  when  the  microwave  frequency  has  swept  beyond  exact 
resonance  by  approximately  the  half  width  of  the  hyperfine  line  shape  observed 
in  slow  passage  (see  Fig.  3).  Thus,  if  we  consider  that  the  transition  is 
power-broadened  so  that  it  has  a  FWHM  of  2u>p  then 

t  -  w  /P  (21) 

m  1 

Furthermore,  this  value  of  t  is  consistent  with  the  standard  physical  picture 

m 

of  ARP.  Consider  ARP  in  the  rotating  coordinate  system  in  the  absence  of 
damping.  As  Fig.  7  shows,  the  atomic  spin  has  a  precessional  frequency  $ : 

?i  -  Bt  z  +  aij  x  (22) 

where  the  projection  of  $  onto  the  z  axis  determines  the  probability  of  find¬ 
ing  the  atom  in  a  specific  spin  state.  The  important  point  to  note,  however, 
is  that  the  projection  of  ft  onto  the  z  axis  only  changes  appreciably  near 
resonance.  When  the  frequency  has  swept  beyond  resonance,  the  population 
reversal  is  nearly  3/4  complete. 

Substituting  Eq.  (21)  into  Eq.  (20)  results  in  an  approximate  expression 
for  the  maximum  degree  of  population  reversal.  To  determine  the  values  of  the 
rapidity  and  adlabatlclty  parameter,  we  need  only  find  the  solution  to 

dlCl<tm>|2 

- — - -  0  for  a  fixed  Rabi  frequency  and  relaxation  rate: 

dp 

*0  *1 +  T"7  )  (23) 

and 

«o  +TT>  <24) 

These  values  are  also  plotted  in  Fig.  5  as  the  dotted  curve,  where  we  have 
Identified  y  with  the  transverse  relaxation  rate. 


In  the  Rotating  Coordinate  System  the  Spin  Vector  Precesses 
about  an  Effective  Field  Hej£,  as  Shown  In  (a);  the  Pre- 
cesslonal  Frequency  Is  ft,  as  Shown  in  (b);  and  (c)  Depicts 
the  Preces6ional  Angle  Swept  Out  during  the  Passage  through 
Resonance,  ti,  as  the  Effective  Field  is  Swept  from  -z  to 


V.  DISCUSSION 


A.  COMPARISON  OF  THE  EXPERIMENTAL  RESULT  AND  THE  FULL  87Rb  DENSITY- 
MATRIX  CALCULATIONS 

From  Fig.  5  we  see  that  the  agreement  between  our  experimental  results 
and  the  ®^Rb  theoretical  model  is  quite  good,  and  this  indicates  that  at  least 
on  a  formal  level  we  have  an  accurate  physical  description  of  the  rapidity  and 
adiabaticity  parameters.  In  particular,  the  theory  supports  our  observation 
of  the  near  constancy  of  the  adiabaticity  parameter,  and  yields  the  numerical 
relation 


w j2/e  ~  1.5  (25) 

for  this  range  of  normalized  Rabi  frequency.  This  result  then  implies  the 
inverse  proportionality  between  the  rapidity  parameter  and  the  normalized  Rabi 
frequency: 


o»iT2 


1.5 

“lT2 


(26) 


B.  COMPARISON  OF  THE  87Rb  DENSITY-MATRIX  CALCULATIONS 

AND  THE  SPIN-1/2  NUMERICAL  SOLUTIONS 

For  o>jT2  >  3  the  numerical  spin-1/2  calculations  and  the  ®^Rb  calcula¬ 
tions  show  excellent  agreement.  In  retrospect  this  should  not  have  been  too 
surprising;  when  the  conditions  for  ARP  are  well  satisfied,  the  extra  Zeeman 
sublevels  of  the  ®^Rb  ground  state  are  essentially  decoupled  from  the  two 
states  interacting  with  the  microwave  field,  and  the  *^Rb  atom  behaves  like  a 
two-level  system.  Thus,  it  is  more  than  coincidental  that  the  agreement 
becomes  good  when  €q  becomes  less  than  unity  (l.e.,  when  the  passage  through 
resonance  becomes  rapid).  However,  the  disagreement  for  <  3  implies 

qualitatively  different  physics  when  the  extra  Zeeman  levels  of  the  ^^Rb 
ground  state  contribute  to  the  ARP  phenomenon.  In  this  range  of  normalized 
Rabi  frequency,  intramultiplet  relaxation  becomes  Important. 


The  agreement  between  the  two  calculations  also  implies  that  the 
constancy  of  the  adiabaticity  parameter  is  a  general  phenomenon  of  ARP.  To 
support  this  claim  we  note  the  recent  experiment  of  Nakanishi  et  al.,^  where 
the  degree  of  population  reversal  in  optical  ARP  was  measured  as  a  function  of 
the  resonance  frequency  sweep  rate.  The  values  of  the  rapidity  and  adiaba¬ 
ticity  parameters  obtained  from  their  experimental  results  are  represented  by 
X's  in  Fig.  5.  Furthermore,  Ernst  and  Anderson^  have  performed  numerical 
calculations  to  determine  the  sweep  rate  that  maximizes  the  absorption  mode 
signal  in  NMR,  and  have  obtained  a  similar  constant  for  the  adiabaticity 
parameter.  Finally,  Antonlewicz  and  Ruijgrok^  have  suggested  that  in  the 
absence  of  longitudinal  relaxation,  but  in  the  presence  of  transverse  relaxa¬ 
tion,  the  adiabatic  condition  for  population  reversal  is  weaker  than  that 
defined  by  Eq.  (2): 

(uJ/B)2  »  1  (27) 

C.  QUALITATIVE  ANALYSIS  OF  THE  RESULTS 

Our  approximate  formulas  for  the  rapidity  and  adiabaticity  parameters 
show  reasonable  agreement  with  the  numerical  calculations,  considering  the 
severity  of  the  approximations  used  in  deriving  them,  and  should  therefore 
contain  qualitative  information  concerning  the  adiabaticity  parameter's  nearly 
constant  value.  Rewriting  Eq.  (20),  we  see  that  the  maximum  degree  of 
population  reversal  is  given  by 

lCl(tm)|2  ~  l1  ”  «XP<“  *h0/2)J  exp(-  £0)  (28) 

which  can  be  thought  of  as  the  product  of  two  probabilities.  The  first  term 
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is  taken  as  the  probability  for  an  adiabatic  passage, and  the  second  term  as 
the  probability  for  relaxation  not  to  occur  during  the  passage: 

P.dl.b.tlc  '  11  '  ’V2)1  <2,a) 


(29  b) 


and 

From  the  conditions  for  ARP  as  stated  in  Eq.  (1),  one  would  expect  the  degree 
of  population  reversal  |Cj(tm)|^  to  approach  unity  for  some  value  of  the  sweep 
rate  as  the  Rabl  frequency  becomes  large,  and  from  Eq.  (28)  this  implies  the 
following  two  limits: 


'  relax 


~  [ 1  -  exp(- 


n(/“lT2^ 


and 


lim 

(0)^2  ♦  *) 


'  adiabatic 


(30a) 


lim 

(i»)jT2  ♦  «) 


relax 


(30  b) 


From  Eq.  (29a),  the  first  of  these  limits  implies  that  nQ  is  an  increas¬ 
ing  function  of  the  normalized  Rabi  frequency: 


lim 


<»1T2 


(31) 


However,  from  Eq.  (29b)  and  the  second  limit  we  have 

lim  n, 


0 


(“lT2*->  U.T: 


0 


r2 


(32) 
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Now,  using  l'Hospital's  rule,  we  can  obtain  the  limiting  value  for  the  rate 
at  which  the  adiabaticlty  parameter  changes  as  a  function  of  normalized  Rabl 
frequency: 


lim  n0  lim  (dn0/d(u1T2)J 

(rn.T.  ♦  •)  (w,T«)  (w.T„  ♦  •)  [d(w,T„)/d(w,T„)] 


lim  driQ 
(w . •)  d(w.T-) 


0  (33) 


Thus,  though  the  adlabaticlty  parameter  is  an  increasing  function  of  the 
normalized  Rabi  frequency,  we  would  expect  it  to  Increase  slowly;  furthermore, 
it  would  appear  that  this  is  due  to  the  necessity  of  minimizing  the  amount  of 
relaxation  that  occurs  during  passage* 

In  order  to  understand  the  physical  implication  of  the  nearly  constant 
value  we  have  obtained  for  the  adlabaticlty  parameter,  we  should  review  the 
physical  content  of  the  ARP  adiabatic  condition.  Referring  to  Fig.  7,  which 
shows  the  motion  of  a  spin-1/2  system  in  the  rotating  coordinate  system,  the 
adiabatic  condition  states  that  the  precessional  frequency  ft  should  be  large 
enough  to  allow  the  spin  to  follow  changes  in  the  direction  of  the  effective 
field.  Put  differently,  the  precessional  angle  swept  out  by  the  spin  in  a 
time  Interval  At  should  be  of  sufficient  magnitude  so  that  the  spin  can  track 
the  corresponding  change  in  the  direction  of  the  effective  field. 

We  now  denote  by  A$  the  angle  through  which  the  spin  precesses  in  a  time 

At  (A$  -  ft  At)  and  consider  a  power-broadened  transition  of  full  width  2  «j. 

The  time  it  takes  to  cross  the  resonance  is  then  At  ■  2  Uj/B,  and  in  that  time 

the  effective  field  has  rotated  by  90® .  Near  resonance  we  can  assume  that  the 

precessional  frequency  is  a  minimum  (ft  »  w.)  and  that  the  precessional  angle 

9  A 

swept  out  by  the  spin  is  A$r  ~  2(^/0).  If  A$r  is  to  be  large  enough  so  that 
the  spin  tracks  the  90°  change  in  the  direction  of  the  effective  field,  then 
the  adlabaticlty  parameter  must  be  larger  than  some  minimum  value. 

From  the  fact  that  the  adlabaticlty  parameter  is  constant  over  such  a 
wide  range  of  Rabi  frequency,  where  we  know  the  spin  is  following  the  effec¬ 
tive  field  (i.e.,  we  can  see  the  effects  of  the  population  reversal),  we  must 
assume  that  the  minimum  value  of  the  adlabaticlty  parameter  is  ~1.5.  This 
would  then  imply  that  the  minimum  precessional  angle  necessary  for  the  spin  to 
track  a  90°  change  in  the  direction  of  the  effective  field  is  three  radians, 
or  approximately  a  half  revolution.  This  is  a  very  surprising  result,  since 
we  would  have  imagined  that  many  revolutions  would  be  required  to  track  a  90® 
change  in  the  direction  of  the  effective  field;  yet  the  above  argument 
indicates  otherwise. 


VI.  SUMMARY 


In  conclusion,  we  have  investigated  ARP  in  the  0-0  hyperfine  transition 
of  *^Rb.  In  particular,  we  have  measured  and  calculated  the  value  of  the 
adiabaticity  parameter  that  maximizes  the  degree  of  population  reversal  as  a 
function  of  normalized  Rabi  frequency.  We  find  that  this  value  is  nearly 
constant  over  the  range  of  normalized  Rabi  frequencies  considered,  with  a 
value  of  ~1.5.  The  reason  for  its  constancy  seems  to  be  due  to  the  necessity 
of  minimizing  the  amount  of  relaxation  that  occurs  during  the  passage. 
Furthermore,  the  result  implies  that  for  a  two-level  system  a  precession  of 
only  half  a  revolution  about  the  effective  field  is  sufficient  for  the  spin 
vector  to  follow  a  90°  change  in  the  direction  of  the  effective  field. 

Through  our  calculations  and  comparisons  of  the  present  results  with  other  ARP 
experiments,  we  have  demonstrated  the  applicability  of  our  conclusions  to  ARP 
in  general. 
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APPENDIX:  THE  RELATION  OF  THE  TRANSMISSION  LINE  SHAPE  LINEWIDTH  TO 
THE  PEAK  AXIAL  MICROWAVE  RABI  FREQUENCY 


In  the  present  experiment  we  used  the  transmission  line  shape  linewidth 
to  measure  the  microwave  Rabi  frequency.  However,  since  the  experiment 
employed  a  TEqjj  cavity  with  its  concomitant  mode  structure,  and  because  it  is 
not  valid  to  assume  that  the  motion  of  the  atoms  performs  a  complete  averaging 
of  the  microwave  field, ^  care  must  be  exercised  in  interpreting  the 
measurement  of  the  Rabi  frequency  by  this  method.  As  stated  in  the  text,  the 
radial  distribution  of  the  microwave  power  can  be  considered  as  essentially 
uniform.  However,  the  longitudinal  dependence  of  the  microwave  field  has  the 
form  sin(^) ,  where  L  is  the  length  of  the  cavity  and  z  is  the  axial  position 
in  the  cavity  measured  from  the  front  face.  In  this  appendix  we  present  an 
argument  to  show  that  in  our  experiments  the  transmission  line  shape  linewidth 
provides  a  reasonable  measurement  of  the  peak  axial  microwave  Rabi  frequency. 


For  simplicity,  consider  the  transmission  line  shape  for  an  optically 
pumped  8pin-l/2  atom.  It  is  quite  easy  to  show  that  the  transmitted  light 
Intensity  as  a  function  of  microwave  detuning  A  is  given  by  an  expression  of 
the  form 


L 

1(A)  -  IQ  exp[-Na  /  c(A,z)  dz] 


(Al) 


where  N  and  o  are  the  number  density  and  absorption  cross  section,  respec¬ 
tively.  ((A,z)  is  the  fraction  of  atoms  in  the  absorbing  state  as  a  function 

O 

of  axial  position  and  detuning: 


C(A,z) 


1  r—l—')  r  i  ♦  * 

2  ''R  +  yJ  1  y 


(R  + 


« 2  2 
y)  +  A 


(A2) 


I 

where  R  and  y  are  the  pumping  and  relaxation  rates,  respectively.  The  signal 
is  the  change  in  the  transmitted  intensity  as  the  microwaves  are  swept  across 
resonance:  S  ■  [I(A  *  »)  -  1(A)]. 
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If  the  vapor  is  optically  thin,  then  the  signal  takes  on  a  fairly  simple 


8  -  I  (rr?)  t - r*— 5 - j)  <**  <A3> 

Furthermore,  by  having  an  optically  thin  vapor  we  can  assume  that  the  pumping 
rate  is  constant  over  the  length  of  the  cell,  which  thus  leaves  the  Rabi  fre¬ 
quency  as  the  only  variable  dependent  on  the  axial  position.  Replacing  u 
by  a>lp  sin(irz/L)  (i.e.,  the  atoms  are  rigidly  frozen  in  place  because  of  the 
presence  of  a  buffer  gas),  we  have 


L  u),  2  sin2  (irz/L) 

S  ~  /  [ - 2~*P-2 - 2 - 2 - ’]  dz 

0  (R  +  y)  +  A  +o).  sin  (rtz/L) 

IP 


Evaluating  the  integral  yields 


S  ~  it  (1  -  /x) 


(R  +  y)2  +  A2 

(R  +  Y>2  +  A2  +  o).  2 
IP 


Though  the  signal  line  shape  expressed  by  Eq.  (A5)  does  not  have  a 
particularly  appealing  form,  it  is  quite  easy  to  show  that  when  conditions  are 
chosen  so  that  x  approaches  unity,  the  signal  line  shape  is  Lorentzlan: 


From  Eq.  (A6)  it  is  clear  that  Eq.  (A7)  is  a  good  approximation  of  the  predic- 

2  2 
ted  signal  when  lu,  /(R+y)]  «  1  and/or  is  in  the  wings  where  (w.  /A)  «  1. 

1  p  1  p 

Furthermore,  it  is  obvious  that  the  linewidth  of  the  approximating  Lorentzian 
is 

aJ/2  -  (HWHM)2  -  (R  +  y)2  +  <*Jlp2  (A8) 

In  Fig.  A-l  we  show  the  line  shape  predicted  from  Eq.  (A5)  under  the  high 
microwave  power  conditions  of  the  present  experiment,  along  with  the  approx¬ 
imating  Lorentzian  of  Eq.  (A7).  Though  the  linewidths  of  the  two  line  shapes 
differ  by  roughly  a  factor  of  two,  it  should  be  noted  that  there  is  good 
agreement  in  the  wings  where  Eq.  (A7)  is  valid.  Therefore,  if  the  above 
discussion  of  the  signal  line  shape  contains  a  complete  description  of  the 
relevant  physics,  then  we  must  conclude  that  the  linewidth  of  the  transmission 
signal  does  not  represent  a  good  measure  of  u^p*  rat^er»  one  should  obtain 
this  quantity  by  examining  the  line  shape  wings.  However,  if  we  consider  more 
closely  the  physical  assumptions  made  in  writing  Eq.  (A2),  we  realize  that  the 
analysis  of  the  line  shape  presented  in  this  appendix  does  not  consider  the 
boundary  conditions  on  c  Imposed  by  the  walls  of  the  absorption  cell.^  It 
might  be  reasonable  to  ignore  these  boundary  conditions  when  R  »  y  (i.e., 
when  the  distribution  of  hyperfine  polarization  is  nearly  constant  over  the 
entire  length  of  the  cell,  abruptly  falling  to  zero  at  the  walls),  but  one 
cannot  expect  this  to  be  a  reasonable  assumption  under  the  conditions  of  the 
present  experiment,  where  R  >  y.  In  fact,  under  the  conditions  of  the  present 
experiment  one  would  expect  the  spatial  regions  furthest  from  the  walls  to 
contain  the  greatest  degree  of  hyperfine  polarization,  so  that  these  regions 
should  be  weighted  more  heavily  in  the  Integral  of  Eq.  (A4). 

Figure  A-2  shows  the  experimental  slow-passage  line  shape  under  condi¬ 
tions  of  high  microwave  power  along  with  a  Lorentzian  fit  of  the  data  in  the 
wings.  (The  data  are  a  signal  average  of  1000  slow  passages  through 
resonance.)  As  can  be  seen,  the  agreement  between  the  experimental  line  shape 
and  the  approximating  Lorentzian  is  much  better  than  Eq.  (A5)  would  lead  one 
to  expect.  As  mentioned  above,  the  reason  for  the  better  agreement  lies  in 
the  fact  that  at  low  light  Intensity  most  of  the  signal  comes  from  the  spatial 
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Fig.  A-2. 


The  Slow-Passage  Experimental  Line  Shape  (solid  curve)  for 
High  Microwave  Power ,  and  a  Least-Squares  Lorentzian  Fit 
(dotted  curve)  to  the  Wings  (hatched  areas).  The 
reasonable  agreement  between  the  two  curves  Indicates  that 
the  HWHM  of  the  experimental  line  shape  is  a  fairly 

accurate  measure  of  u,  . 

IP 
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regions  furthest  from  the  walls.  It  is  a  coincidence  of  the  TEq^j  mode 
structure  that  the  regions  expected  to  contain  the  smallest  values  of  5 
correspond  to  the  regions  where  k)^(r,z)  =  “ip*  Thus,  we  make  the 
approximation  that  the  HWHM  of  the  experimental  line  shape  is  indicative 


of  oi,  ,  and  Fig.  A- 2  indicates  that  this  will  lead  to  <  10%  error,  which  is 
ip  ~ 

insignificant  on  the  log  plot  of  Fig.  5. 


LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experlaental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  adentlflc  advances  to  new  military  space  ays  teas.  Versatility 
and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  person¬ 
nel  in  dealing  with  the  aany  problems  encountered  in  the  nation's  rapidly 
developing  apace  systeaa.  Expertise  In  the  lstest  scientific  developeents  Is 
vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The  labora¬ 
tories  that  contribute  to  this  research  are: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chealstry,  environmental  hazards,  trace  detection;  spacecraft  structural 
mechanics,  contamination,  thermal  and  structural  control;  high  temperature 
thermoaechanlcs,  gas  kinetics  and  radiation;  cw  and  pulsed  laser  development 
Including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control, 
atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency 
standards,  and  environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  space borne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence  .and 
microelectronics  applications. 

Electronics  Research  Laboratory;  Microelectronics,  GaAs  low  noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro¬ 
optics;  communication  sciences,  applied  electronics,  semiconductor  crystal  and 
device  physics,  radiometric  Imaging;  millimeter  wave,  microwave  technology, 
and  RF  systems  research. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal 
matrix  composites,  polymers,  and  new  forms  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  reliability;  fracture  mechanics  and 
stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures  as  well  as  In  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnatospherlc,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  Interactions,  magnetospherlc  plasma  waves;  atmospheric  and 
Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote 
sensing  using  atnospherlc  radiation;  solar  physics.  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
Instrumentation. 


